Theoretical investigation of TbNi 5 A Cu A optical properties 
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In this paper we present theoretical investigation of optical conductivity for intermetallic TbNij_ v Cu v series. 
In the frame of LSDA+U calculations electronic structure for x = 0, 1, 2 and on top of that optical conductivities 
were calculated. Disorder effects of Ni for Cu substitution on a level of LSDA+U densities of states (DOS) 
were taken into account via averaging over all possible Cu ion positions for given doping level x. Gradual 
suppression and loosing of structure of optical conductivity at 2 eV together with simultaneous intensity growth 
at 4 eV correspond to increase of Cu and decrease of Ni content. As reported before [Knyazev et ah, Optics and 
Spectroscopy 104, 360 (2008)] plasma frequency has non monotonic doping behavior with maximum at x = 1. 
This behavior is explained as competition between lowering of total density of states on the Fermi level N(E F ) 
and growing of number of carriers. Our theoretical results agree well with variety of recent experiments. 

PACS numbers: 



I. INTRODUCTION 

Intermetallic compounds of RNis-type have been attract- 
ing considerable interest initiated by their remarkable physic - 
ochemical properties and potential applications as materials 
with high hydrogen storage capacity-^. Their physical prop- 
erties, which are characterized by a large variety of magnetic 
structures and electronic peculiarities, have been actively in- 
vestigated, both theoretically and experimentally. Substitution 
of one rare-earth metal R with another one modifies signifi- 
cantly these characteristics because of crystal field effects and 
changes of exchange interaction value between the conduction 
and 4f-localized electrons. 

On the the other hand substitution of Ni for some other tran- 
sition metal also leads to interesting physical phenomena. For 
example, series of isostructural RNi5_ A Cu A alloys was found 
experimentally to have non monotonic concentration depen- 
dence of the magnetic susceptibility, Curie temperature Tc, 
electronic specific heat and resistivity 3-7 . A number of studies 
indicate a direct correlation between this anomalous behavior 
of the parameters and band structure evolution with copper 
content increase 4,5,8 .. To explain these experimental data one 
needs more detailed investigation of electronic structure of the 
RNis A Cu A series for different x. 

There are several experimental and theoretical works on 
the RNiCu compounds. X-ray photoemission study of 
LaNi5_ v Cu A solid solution shows presence of Cu-3d band, 
which is located almost 2 eV below the 3d band of Ni and 
is weakly hybridized with the latter one 8 . The calculation of 
Ref. |U predicts the existence of a wide (~ 1,5 eV) peak in the 
electron density of states of the YNi5_ A Cu A alloys at energies 
2.5-4 eV below Fermi level E^ , which is related to the 3d elec- 
trons of Cu. The evolution of frequency dependencies of the 
optical conductivity of the TbNi5_ A Cu A compounds upon sub- 
stitution of Ni for Cu atoms investigated in Ref. LKJ. There it 
was observed that a new structure arises in the optical spec- 
tra (a broad absorption maximum peaking at 4 eV) while the 
Cu content increases. Detailed investigation of crystal struc- 
ture and thermodynamic properties from both experimental 
and band structure calculation sides is done in Ref. 9 (except 



for the optical properties). 

Within this study we continue investigation of the electronic 
structure of TbNi5_ A Cu A isostructural alloy series with x — 0, 
1, 2 by carrying out self-consistent ab initio LSDA+U calcu- 
lations. Disorder effects are taken into account in a combi- 
natorial way with respect to possible number of Cu positions 
for given doping. Corresponding averaged over the Cu posi- 
tions spin resolved densities of states are calculated. Follow- 
ing it interband contribution to the theoretical optical conduc- 
tivity is determined. Drude (intraband) contribution is taken 
as a Lorentzian with experimentally determined parameters. 
Comparison of our theoretical results with recent optical ex- 
periment is presented. 



II. CALCULATION OF THE ELECTRONIC STRUCTURE 

The TbNis compound crystallizes into a hexagonal struc- 
ture of the CaCus type with a space group P6/mmm and lat- 
tice parameters for x=0 — a=b=4.8988 A and c=3.96 A (see 
Fig. [Qh x=l — a=£=4.9126 A and c=3.9901 A and x=3 
_ a =b=A.935\ A and c=4.0113 A. 9 The nickel has two in- 
equivalent crystallographic positions Nil(2c) (1/3,2/3,0) and 
Ni2(3g) (1/2,0,1/2), Tb has a position 1(a) (0,0,0). In Fig.Q] 
Nil ions are represented by violet circles and located in the 
same "layer" as Tb ions (green circles, Fig. [TJ. The Ni2 ions 
(red circles, Fig.Q]) also form a "layer". 

A calculations of the electronic structure of TbNi5_ t Cu A 
(x = 0, 1, 2) were performed within the linear muffin-tin or- 
bitals (LMTO) method^. The program package TB-LMTO- 
ASA (tight-binding LMTO together with atomic spheres ap- 
proximation ) v. 47 was employed 12 . The subdivision of the 
first Brillouen zone in the reciprocal space is (10 10 10) with 
132 irreducible k-points. The atomic sphere radius for Tb 
was obtained to be about 3.47 a.u., for nickel and copper ions 
about 2.63 a.u. The orbital basis consists of 6s, 5d and 4f 
muffin-tin orbitals for Tb, and 4s, 4p and 3d for Cu and Ni. 
LSDA+U calculations on top of that were done as proposed in 
Ref.[l4 with Coulomb interaction parameters for Tb(4f) states 
U — 5.4 eV and J = 0.7 eV4^ Also recently we showed that 
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FIG. 1: (Color online) TbNis crystal structure. Tb ions 
cles, Nil,Ni2 - violet and red circles correspondingly. 



green cir- 



Coulomb interaction for Fe ions in intermetallides is almost 
completely accounted within LSDA by itself. 13 

To introduce Cu into the TbNis system one or two Ni ions 
were substituted with Cu. To mimic disordered alloy we per- 
form calculations for different possible positions of copper 
ions for chosen doping level x. Assuming that each pos- 
sible Cu position has the same weight we then average ob- 
tained partial densities of Ni and Cu-3d states with following 
ratios for corresponding Ni Wickhoff positions: for x — 1 
i.e. one copper ion 2(2c):3(3g) and x = 2, two Cu ions, 
l(2c):3(3g):6(2c+3g). 

In our LSDA+U calculations collinear ferromagnetic or- 
der of the local magnetic moments on all lattice sites was 
chosen. Moreover the terbium and nickel sublattices are 
taken to be ferromagnetic with respect to each other. In this 
compound main magnetic moment is localized on terbium 
ions while nickel ions have a small magnetization because 
of almost completely occupied Ni-3d shell in accordance 
with Ref. |2j. For parent TbNis system local magnetic mo- 
ments in units of Bohr magneton are following: on Tb(5.8), 
Nil(0.27) and Ni2(0.28). Total magnetic moment of TbNi 5 is 
7.21 Bohr magneton. With Ni for Cu substitution in TbNis 
depending on Cu position total magnetic moment goes down 
about 10-15%. These results agree well with experimental 
data of Ref. 9. Local magnetic moments on nickel ions be- 
come two to three times smaller. Local magnetic moment on 
terbium insignificantly grows, while copper stays practically 
non magnetic. 

The Fig. [2] shows the total (upper panel) and partial Ni(3d) 
(middle panel) and Cu(3d) (lower panel) DOSes for up (T) and 
down (X) projections of local spin moments for the ferromag- 
netic TbNi5_ A -Cu v (x=0,l,2) compounds. Structure of these 
DOS for different x is rather similar to each other and agrees 
with Ref. [9] for stoichiometric TbNis. It is seen that the main 
spectral weight is located below the Fermi level. The 3d states 
of Nil and Ni2 form the wide bands extending from 4.0 to 0.5 
eV. The 3d states of Cul and Cu2 bands for the compounds 
with x=l,2 are located in the energy range -4 — 6 eV below 
Ef . Their width is considerably smaller with respect to Ni-3d 
states. The 4f-states of Tb ion for majority spin are located 
below -6 eV and for minority at about 1.5 eV. 



III. OPTICAL CONDUCTIVITY RESULTS AND 
DISCUSSION 



Calculated LSDA+U band structures presented in the pre- 
vious section are used to interpret experimental optical data 10 
for the systems under consideration. In order to calculate 
theoretical interband optical conductivity o,i, eor we applied a 
rather simplified technique 16 using approximation, that the di- 
rect and indirect (involving phonons) interband transitions are 
equally probable. Namely we computed o,i ieor like the inte- 
gral function based on the convolution of the density of states 
located both below and above the Ef. The electronic transi- 
tions of the d— >f type are not observed in the calculation due to 
their low probability. This technique was already successively 
applied by us to intermetallic R2Fen compounds^. 

From total LSDA+U DOS presented on upper panel of 
Fig. [2] one can observe decrease of total DOS on the Fermi 
level N(Ef) with Cu doping. It happens because for x=0 
N(Ef) consists mostly of Ni-3d states. With doping number 
of Ni ions goes down while newly introduced Cu ions do not 
contribute to the N(Ep) (see lower panel of Fig. [2j». 

Figure|3]displays corresponding theoretical optical conduc- 
tivity a,heor{oS) (solid lines) together with experimental fre- 
qu< 
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FIG. 2: (Color online) LSDA+U calculated DOS for TbNi 5 -. v Cu v . 
Upper panel - total DOS; middle panel - Ni-3d states; lower panel - 
Cu-3d states. In all graphs solid black line represents TbNis; dashed 
blue lines - averaged DOS for TbNi 4 Cui; dash-dot red line - aver- 
aged DOS for TbNi3Cu2. Fermi level corresponds to zero. 
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FIG. 3: (Color online) Comparison of experimental optical conduc- 
tivity (dots) and theoretical one (solid line in arbitrary units) for dif- 
ferent doping levels x=0, 1,2 (from top to bottom) for TbNi5_ x Cu x . In 
the inset of upper panel Drude peak and corresponding experimental 
parameters of the Lorentzian are shown. 

for TbNig_ t CUj. (x = 0, 1, 2 from top to bottom). In the energy 
region below ~ 0.5 eV the behavior of o(w) for all composi- 
tions sharp Drude peak is determined experimentally. Corre- 
sponding height and width of Drude peak (127xl0 14 s and 
0.280 eV for x = 0, 156xl0 14 s" 1 and 0.278 eV for x = 1, 
130xl0 14 s" 1 and 0.288 eV for x = 2 ) were extracted from 
these curves and then used to draw theoretical one (see inset 
of Fig. 0). As follows from above mentioned static conductiv- 
ity (Drude peak height) it goes down with doping. It agrees 
with lowering of N(E F ) in our LSDA+U calculations since 
conductivity is proportional to N(Ef). In its turn plasma fre- 
quency is also proportional to N(Ep) and thus should decrease 
with Cu doping. On the other hand substitution of Ni for Cu 
adds one electron to the system and increases number of car- 
riers. Plasma frequency squared is proportional to the number 
of carriers (in the simplest Drude-like picture) and then should 
grow with Cu doping. These two competing tendencies lead 
to the non monotonic doping behavior of plasma frequency 
recently observed experimentally in Ref. |T(J 

At energies above 0.7 eV the shape of o(w) dispersion in- 



dicates the dominant role of interband absorption. The wide 
absorption region is characterized by structures, which inten- 
sities and localizations depend on the compound composi- 
tion. The optical conductivity of the TbNis binary alloy (upper 
curve) exhibit two peaks at photon energies 1.15 and 1.9 eV 
respectively. In the TbNiztCu compound the o(oj) spectrum is 
almost the same; we can only note that the minimum between 
these structures is less pronounced. With a further increase 
in the Cu content (x = 2) a new fairly broad absorption band 
arises at energies 3-4.5 eV. 

Corresponding theoretical optical conductivity curves o r / leor 
shows evolution upon substitution of nickel with copper sim- 
ilar to experimental one. Theoretical calculations reproduce 
qualitatively the basic features of the experimental spectra, 
namely, two-peak structure at 1-2.5 eV, associated with the 
Ni 3d interband transitions. The depth about 1.5 eV is getting 
closed. High energy tail above 3 eV becomes smoother with 
doping increase and for x — 2 one can see the formation of the 
"Cu-states" peak at about 5 eV. Somewhat higher energy posi- 
tion of leading interband absorption edge (about 1 eV) might 
come from uncertainties in computation of U value. 



IV. CONCLUSION 

In this work we perform ab initio LSDA+U computa- 
tions of electronic structure of intermetallic isostructural al- 
loys TbNi 5 _ A .Cu, (x =0, 1, 2). Disorder effects of Ni for 
Cu substitution were accounted in a combinatorial way. Cor- 
responding averaged over all possible Cu positions for given 
doping level x DOS are further used to compute theoretical 
optical conductivity. Here we is theoretically explain how ex- 
perimental optical conductivity line shape changes with dop- 
ing. Analyzing LSDA+U calculated DOS one can find that 
responsible for the structures of optical conductivities at 1-2 
eV electronic states are predominantly Ni-3d states, while at 
~4-5 eV there are Cu-3d states. Amount of latter ones in- 
creases upon doping giving rise to the structure at 4-5 eV in 
optical conductivity. One should note also good semiquanti- 
tative agreement between theoretical and experimental data. 
Finally in this work we give explanation of the experimen- 
tally observed maximum of plasma frequency for x=l— as 
competition of lowering of N(Ep) and simultaneous growing 
of number of carriers with Ni for Cu substitution. 
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